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Abstract 

In this paper we present an equilibrium value based framework for solving SDPs via the multiplica- 
tive weight update method which is different from the one in Kale's thesis |Kal07]. One of the main 
advantages of the new framework is that we can guarantee the convertibility from approximate to exact 
feasibility in a much more general class of SDPs than previous result. Another advantage is the design of 
the oracle which is necessary for applying the multiplicative weight update method is much simplified in 
general cases. This leads to an alternative and easier solutions to the SDPs used in the previous results 
QIP(2)CPSPACE |JUW09| and QMAM=PSPACE |JJUW09| . Furthermore, we provide a generic form of 
SDPs which can be solved in the similar way. By parallelizing every step in our solution, we are able to 
solve a class of SDPs in NC. Although our motivation is from quantum computing, our result will also 
apply directly to any SDP which satisfies our conditions. 

In addition to the new framework for solving SDPs, we also provide a novel framework which im- 
proves the range of equilibrium value problems that can be solved via the multiplicative weight update 
method. Before this work we are only able to calculate the equilibrium value where one of the two convex 
sets needs to be the set of density operators. Our work demonstrates that in the case when one set is the set 
of density operators with further linear constraints, we are still able to approximate the equilibrium value 
to high precision via the multiplicative weight update method. 



The work was completed when the author was visiting the Institute for Quantum Computing , University of Waterloo. 
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1 Introduction 



Semidefinite programming (SDP) is a relatively new field of optimization which grew up in the 1990s 1AU951 
BV04. Lov03 VB96. dK02[. Despite of the short time since its introduction, SDP has proved useful in many 
different contexts. Especially there are many applications of SDPs in the field of theoretical computer sci- 
ence, like the design of approximation algorithm BVazOll and the recent application to Unique Game Con- 
jecture JSiinyRagOS) and simulating quantum complexity classes IGut05l IGW051 IGW071 [TTUW09J [TUW09J 
IKW00[IWat09bT 

SDPs are in fact a special case of conic programming and there exists polynomial algorithms to solve 
any SDP instance (like the interior point method |Ali951ldK02ll ). Thus, any problem which can be mod- 
eled or approximated as a SDP is considered to have efficient solutions. However, on the other side, our 
understanding of SDPs is far less than our understanding of Linear programming(LP), another typical op- 
timization method used a lot in practice. One thing is about the running time : although any SDP instance 
can be solved in polynomial time theoretically, it is much slower than LP's solution in practice. Moreover, 
the generic algorithm for SDP are always used as a black box. Hence it is rarely seen that we can employ 
the duality or the structure of SDPs while the duality of LPs inspires lots of new algorithm designs. 

A generic primal-dual method for SDP problem (under certain conditions) was introduced by Arora 
et.al. |AHK05al]AHK05bllAK07 | to overcome the difficulties mentioned above. The generic method exploits 
a generic framework (or meta-algorithm) called the multiplicative weight update method. Similar frameworks 
were studied flFleOOl |FS99l IGK981 IPST911 |Kha04[ IWK061 |You95H for different purposes before . This new 
generic method turns out to be very useful and successful in several contexts. In the paper |AK07] (see 
also Kale's Phd thesis |Kal07|) where it was originally proposed, this generic method improves the upper 
bounds of running time of many approximation algorithms. In addition to that, this generic method was in- 
troduced by Watrous et. at to the field of quantum computation and successfully proved QIP(2)CPSPACE [ JUW09 [ 
and QIP=PSPACE |JJUW09| . 

A semidefinite program over X and y (shown below) is specified by a triple (Y, A, B) where Y : L [X ) — > 
L (y) is a Hermiticity preserving super-operator and A £ Herm [X) and B 6 Herm (y). This form which 
was obtained in [Wat08J is somehow different from but equivalent to the standard form. Let a, 6 be the 
optimum values of the primal and dual programs respectively. One important property of the semidefinite 
program, called the duality, implies that a < 6 and the equality will hold in some situation. 

Primal problem Dual problem 

maximize: (A,X) minimize: (B, Y) 

subject to: Y(X) < B, subject to: Y*(Y) > A, 
XeFos(X). YePos(J). 

In practice we usually consider the SDPs whose optimum values are within a small range. For those 
SDPs, instead of directly calculating the optimum value, we usually consider the feasibility problem first. 
Once the feasibility problem is solved, we could use the binary search in that range to find the optimum 
value. For any instance of SDP and a guess value c , the feasibility problem is defined to be 

Feasibility Problem 
ask whether: {A, X) > c 
subject to: Y(X) < B, 

x e Pos (X) . 

Intuitively, the primal-dual method in Kale's Thesis [Kal07| for the feasibility problem contains the fol- 
lowing three ingredients. First, update via the multiplicative weight update method. Second, the existence 
of an efficient width-bounded oracle 0\ . Finally, this method will generate an approximately feasible solution 
Y to the dual problem such that (B, Y) < c . By approximate feasibility we mean 

VpeD(X) ,(Y*(Y)-A,p) > -e (1) 
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where the D (X) denotes the set of density operators over the space X and e is some small constant. In 
order to solve the feasibility problem, we need to get an exact dual feasible solution Y such that ( B, Y ) < 
(1 + e)c from this. We will refer this as the convertibility from approximate to exact feasibility. After executing 
a combination of those three ingredients, this method will either return a feasible solution X to the primal 
problem with object function value at least c or a feasible solution Y such that (B, Y) < (1 + e)c. The latter 
case will imply a. < & < (1 + e)c by the duality of SDPs. A detailed description of this procedure can be 
found in Appendix [A3] 

Another important value which can be calculated via the multiplicative weight update method is the 
equilibrium value of zero-sum games [vN28[ and its generalizations (like, ||Haz06l ). Particulary, we consider 
the value A, 

A = minmax fix, y) = maxmin /(x,y) 

for some convex-concave function / (see definition in Appendix |A.3t over X x Y where X, Y are convex 
compact sets. Again, this method involves the update via the multiplicative weight update method and an 
efficient width-bounded but functionally different oracle O2 as main ingredients. However, for equilibrium 
value, there is no requirement for the convertibility from approximate to exact feasibility. Under several 
other conditions as well (see details in Theorem [2] and Appendix \A3\ , such value A can be approximated 
to high precision efficiently. It should be noted here the multiplicative weight update method plays a quite 
different role from the one in the solution to SDPs above. This difference has led to an alternative and easier 
proof of QIP=PSPACE IWulOI . It was also applied to the proof of QRG(1)CPSPACE |JW09| . 

One important advantage of the solutions to SDPs or equilibrium value based on the multiplicative 
weight update method is that we can easily implement the algorithm in parallel. Precisely, this is because 
the fundamental operations of matrices and singular value decomposition of matrices [Gat93| can be im- 
plemented with high accuracy in NC. This trick was widely exploited in the recent progress of quantum 
complexity theory |Wulfj |jJU\^|JUW09||JW09l . 

In this paper we will demonstrate how the equilibrium value can be related to the feasibility problem 
(then, SDPs). It will then suffice to make use of the solution to the equilibrium value to solve the feasibility 
problem. Following this idea, we will provide an alternative and easier solution to the SDPs used for 
QIP(2) [JUW09[ and QMAM |JJUW09] . Moreover, we provide a generic form and conditions under which 
any feasibility problem can be solved in the same way. By parallelizing each step in the solution, we are able 
to show any feasibility problem can be solved in NC as well. Precisely, we consider the feasibility problem 
as follows. 

Feasibility Problem 
ask whether: (A, X) > c 
subject to: Y(X) < B, 

XeD(X). 

This feasibility problem is very similar to the most general version above. The only difference is we 
replace the condition X £ Pos ( X) by X £ D (X). This change corresponds to the trace bound Tr X = R for 
some const R which commonly appears in applications of SDPs. Under certain conditions our constraint 
on X £ D (X) is equivalent to the trace bound Tr X < R. This implies the feasibility problem in our consid- 
eration is very general and could cover many instances in practical use. Although our original motivation 
is from quantum computation, our result also works well for any SDPs which can be converted to the form 
in our consideration. 

The concept of equilibrium values will be related if we imagine a two-player game to solve the feasibility 
problem. Assume there is a primal player who wants to provide you a feasible solution X to prove the 
original problem is feasible. On the contrary, the dual player who wants to disprove the feasibility will try 
to find where the constraints on X are violated. This is different from Kale's method where the disproof of 
the feasibility is by getting some feasible solution to the dual problem and then making use of the duality 
of SDPs. Precisely, we define the following convex-concave function f to capture the two-player game. 



3 



Framework 1. Let function f be 

/(X ,n) = (( c -< A - x > TW _„).n) (2) 

over the set D (X) x T where T = {n : < n < lyec}- the equilibrium value A* be 

A* = min max f (X,n) = max min f(X, FT) 
xeD(X) neT neTxeD(X) 

The relation between the equilibrium value A* and the feasibility of the original problem is captured by 
the following theorem. 

Theorem 1. TTze original problem is feasible if and only if A* < 0. 

This formulation is inspired by similar formulations [Haz06[ used for convex optimization and Theo- 
rem Q] follows easily from the argument in |Haz06| . Given the fact that equilibrium value A* can only be 
calculated approximately, we still need to do the conversion from approximately to exactly feasible solu- 
tions. Due to this reason, we make an important change to the old formulation in [Haz06[. Namely, we 
choose rtobe{]T:0<n< ly^c} rather than the set of density operators over the space y © C. As- 
sume we approximate A* to precision e and the return value implies A* is in an interval containing 0. Or 
equivalently, we encounter the situation 

v n e x,(( c -< A ' x > Y(x) _ B ),n)<e (3) 

and we want to convert X into some exact feasible solution without changing the object function a lot. This 
is very similar to the situation captured by Equation (TJ. However, the difficulty in making this conversion 
happen is different. By applying the approximate feasible result in Equation (U, we can only expect to solve 
the problems |JUW09| |JJUW09| when A = \ x efficiently in general. This is because the density operator 
p in Equation JT) implies the approximate feasibility only holds in term of £oo norm. On the other side, in 
order to remain the object function which is a inner product almost unchanged, the approximate feasibility 
needs to hold in a stronger sense like C\ norm. Otherwise only trivial cases like A = \% can be solved. 

Our new approximate feasibility result in Equation (3J overcomes such difficulty and provides a method 
to do the conversion from the approximate to exact feasibility in much more general cases. Especially we 
will demonstrate when the super-operator Y is partial trace or its generalizations, our approximate feasi- 
bility result works very well to make the conversion happen. This kind of constraints is very powerful 
because it is the only type of constraints we need in many SDPs in quantum computation. In addition to 
that, we still have the freedom to choose T in order to meet the requirement of new types of constraints. 
Such freedom is a huge advantage over the primal-dual method since Equation (T) is the only result can be 
expected from the primal-dual method. 

Another advantage of the Framework [1] is that the oracle O2 which will be required to compute the 
equilibrium value A* can be easily designed. Precisely, since T={n:0<n< ly^c}/ we can simply 
get the other part's spectral decomposition and let T be the projection onto the positive eigenspace of it. 
Finally, under the conditions of Theorem 12 we can implement the whole algorithm in NC. 

Besides the application to the feasibility problem, the idea of Framework [Jean also improve our ability 
on calculating equilibrium over some non-density operator set. Consider some convex-concave function h 
over X x Y where X is the set of density operators but with some constraint O. Precisely, 

Framework 2. If we consider the equilibrium value u, 

u = min max /z(x,y) = maxmin/z(x,y) 

xex yeY yeY xex 

where X = {p : p £ D (X) ,<&(p) < B} and Y is any other compact convex set. It will be useful to consider the 
equilibrium value u' as follows 

u' = min max h{x,y) + a (<&(p) — B, n) = max mmh(x,y) + a (<&(p) — B,H) 

xex' {y,n}eY' {y,n}eY' xex' 

where X = D (X) , Y' = Y x T (T = {IT : < Fl < 1}) and a is any factor. 
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It is easy to see that the idea of Frame work [2] or especially the term (<t>(p) — B, TV) is to penalize when p 
does not satisfy the constraint ^(p) < B. Furthermore, the penalization is weighted according to the factor 
a. By using the game value of one restricted model of one-round quantum refereed game as an example, 
we can demonstrate when there are two promise values of p with large gap, such gap can be transferred 
to the new value p'. Thus, it will be sufficient to calculate the value of p' in order to distinguish between 
two promises of the original value p. This improves the range of problems which can be solved by the 
multiplicative weight update method since so far we can only calculate the equilibrium value when X is 
the set of density operators (up to a factor). This will also give a binary search method for calculating the 
equilibrium value to some precision by artificially assuming two promise values. 

The rest of the paper is organized as follows. Most of the preliminaries can be found in Appendix[A] We 
also leave the lemmas and theorems which will be directly used in Section [2] The two examples for Frame- 
work[T]will be demonstrated in Section|3](QIP(2)) and in Appendix[B](QMAM) respectively. One restricted 
version of one-round quantum refereed game will be discussed in Section [4] to demonstrate the power of 
Framework[2] We will conclude the paper with further discussions and open problems in Section[5] 

There are two points to make clear before the readers move on to the next section. First, when we 
consider the feasibility problem or equilibrium value directly, we always refer the size of the SDP or the 
function with equilibrium value as the input size. However, when we consider the quantum complexity 
classes, the SDP or the function with equilibrium value will have exponential size in term of its actual input 
size \x\. Second, we will not take care of the precision issues with the NC implementation in the main part 
of this paper. Instead, we will assume such implementation can be made exactly and deal with them in 
Appendix ICl 

2 Preliminaries 

In order to make this paper self-contained, we try to provide brief surveys on each topic related to our paper. 
However, most of them will be put in the appendix due to limited space. Precisely, we will introduce the 
fundamentals of quantum information in Appendix lA.il Useful facts on NC and parallel matrix operations 
are stated in Appendix lA.2l The multiplicative weight update method and its application to calculating the 
equilibrium value and SDPs (the primal-dual method) are surveyed in Appendix lA.3l 

2.1 Useful Lemma and Facts in Quantum Information 

The following are some important lemmas about purification and fidelity which are useful in our proof 
later. The proofs of following results are put in Appendix lA.il 

Lemma 1. Given any two density operators p\, pi over the space A, and another density operator <T\ over the space 
A® B such that Trg o~\ = p\, then there exists another density operator cr*i over the space A ® B for which that 
Trg o~2 = p 2 andF(px,p 2 ) =T{(T Xl cr 1 ). 

Lemma 2. In addition to the result in Lemma\l\ we can compute the classical representation ofo~2 as required above 
given the classical representations of p\, p 2 and o~\ in NC where the input size refers to the size of the matrices. 

Lemma 3. Given two density operators p\,p2 € D (A), and their purifications <T\,o~2 6 T)(A<g> B) in space A®B 
respectively , for which T{p\,p2) = ^{o'lrO'i), let 

1 1 

s = — /°2 1| l andt = -\\ai-a 2 \\i 

Then we have the following inequalities 

1 - s < \T\ - t 2 andl-t < Vl ~ s 2 

Lemma 4. Given two density operators p\,p2 € T)(A®B) where p\ represents a pure state, there exists an admiss- 
able quantum channel <T> : L (A) — > L (A) such that <T> <g) 1l(b)(/0i) = Pi if and only ifTr^(pi) = Tr^^)- 
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1. Let £ = £ and T = 16r ^ nD . Also let W' 1 ) =t x ,D = dim (A"). 

2. Repeat for each t = 1, . . . ,T: 

(a) Let jo( f ' = W( f )/Tr W( f ) and let nw be the projection onto the positive eigenspace of S(p^). 

(b) Let M (f) = (N(n( f )) + r\x)l1r and update the weight matrix as follows: 

t 

W( f+1 ) = exp(-£ £ M (T) ) 

T=l 

3. Choose (0 = j= Et=i P an d l et H be the projection onto the positive eigenspace of S(p). Return (p, TV) 
as the approximate equilibrium point and (S(p), TV) as the approximate equilibrium value. 



Figure 1: An algorithm that computes the approximate value and point to precision S. 

Lemma 5. Given any Hermitian operator A such that Tr A = 0, then we have 

1 1 
max (A, VI) = -II All i and min (A.TV) = --\\A\U 
n : o<n<i x ' 2" 11 n:0<n<i N ' 2" " 

2.2 Multiplicative Weights Update Method 

The detailed discussion of the multiplicative weight update method is provided in Appendix lA.3l However, 
we will demonstrate the particular algorithm to calculate the equilibrium value in the form of Equation 
in this section. Precisely, we consider the 

A*= min max/f C ~ ^'^ wf . n \n\ (4) 

The existence of the equilibrium value is implied before. Thus we only need to see when a NC algorithm 
will exist to calculate the A approximately to high precision. To ease the description of the algorithm, let 
S(p) defined to be 

' c-(A,p) 

Y(p) - B 



S(p) 



for any p G T)(X). Choose N(TT) to be the raw loss matrix such that for any p £ T)(X) and IT = 
^ p ^ £ T, we have 

(S( P ) 7 n> = ^N(n)> (5) 

It is easy to see that we can choose N(TT) to be 

N{TT) = - P A + Y*{P) + {pc-{B,P))t x (6) 

Theorem 2. Let the input size, denoted by \x\, be the size of any instance in Equation @. If for any IT £ T the 
N(TT) defined above satisfies || JV(TI) ||oo < r where r = O(polylog(\x\)) and Y(p) can be calculated in NC for 
any p £ D (X), then by using the algorithm in Figure^ we can approximate the equilibrium value to precision 
& = D,(polylog{\x\)) in NC. 

We will leave the proof of Theorem|2]in Appendix IA.3I 
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3 QIP(2) Case 



Now it is our turn to consider a real instance of semidefinite program and apply our framework to solve 
it. Our first candidate is the quantum interactive proof system with two-messages. In this system, after the 
input x is given, the polynomial-time bounded quantum verifier will send one quantum message to an all 
powerful quantum prover and get another quantum message back. Then the verifier will decide whether to 
accept or to reject based on the message sent back from the prover and the qubits kept at his side. The only 
constraint on the all powerful quantum prover is that the prover must operate an admissable quantum 
operation on the quantum message sent to him. The complexity class QIP(2) denotes all the languages 
which can be recognized by the procedure above. Precisely, we have 

Definition 1. Any language L is inside QIP(2) if and only if 

• If x G L, there exists a prover such that the verifier will accept with probability at least c{ \x\ ). 

• If 'x i L,for any prover the verifier will accept with probability at most s(|x|). 
where c(\x\) — s(|x|) = Ci(l/poly{\x\)). 

It is known that QIP(2)CPSPACE |JUW09| by following Kale's way HKal07| to solve SDPs. By contrast, 
we will demonstrate here how our Framework[l]can be applied to this problem. Namely, we provide an al- 
ternative proof of the result QIP(2)CPSPACE. The main difference between our approach and the previous 
approach is that we formulate the problem using the density operators instead of quantum channels and 
we solve the semidefinite program using the new framework. 

Let A4 denote the message's space between the prover and the verifier and V denote the verifier's 
private space. Let us assume the input x is fixed for the following discussion. Without lost of generality, let 
the pure state p\ 6 D (Ai ® V) be the initial state for the input x, namely the state that the verifier prepares 
given input x. The prover will then operate an admissable quantum channel <J> : L (A4) — > L (A4) on part 
of the state p\ and it will result another state pi = <I> <S> 1l(v) The verifier will performance a POVM 
measurement on p 2 to decide whether to accept or to reject. Let R be the POVM which corresponds to 
the case where the verifier accepts. In order to decide whether x £ L, it suffices to solve the optimization 
problem 

max(R,p 2 } s.t. p 2 = <&® l L (V)(Pl) 

where the optimum value is the maximum probability that the verifier accepts given the input x. 

Because of LemmaUl we have p\ and pi are connected by an admissable quantum operation if and only 
if Tr^vj pi = Tr^v) pi. Thus the above optimization problem is equivalent to the following SDP, denoted by 
SDP (I). 

SDP Problem Feasibility Problem 

maximize: {R,P2) ask whether: (R,pi) > c 

subject to: Tr^ (p 2 ) < Tr M (pi), subject to: Tr^ (p 2 ) < Tr M (p{), 

p 2 e D (M ® V) . p 2 ED(M®V). 

3.1 Solution to the Feasibility Problem 

Following the Framework[TJ we consider the feasibility problem above. Precisely, we define 

/,(P,n) = (('-<M TtM(p) _ TtM(fl) ),n) P) 

where p 6 T\ = D (M <g) V) and lie T 2 = {IT : < IT < 1m®c}- Let be the equilibrium value of 
function f\, namely, 

A? = min max fi (p, Tl) = max min fi (p, Tl) 

peTjlleTy r TLeT 2 peTi 
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Base on Theorem [TJ the value of A^ will imply whether the original problem is feasible. In addition to 
that, we will demonstrate how to convert any approximately feasible solution to exactly feasible solution 
without changing the value of the object function a lot. 

Lemma 6. Assume we can calculate the approximate equilibrium value and point of the function f\ to the precision 
S. Let the \\ and {p, IT} be the approximate equilibrium value and point returned by algorithm in Figure^ A* be the 
actual equilibrium value. Then we have 

• if X\> 5, then the original problem is infeasible. 

• Z /A~i < 5, then there exists a feasible solution p such that (R,p) > c — \/2S — S 2 . 

Proof. • If Ai > 5, namely, A^ > Ai — 5 > 0, then due to Theorem [TJ the original problem is feasible. 

• Otherwise, due to Lemma [5j we have 

h(p,P) = max{c - (R,p) ,0} + - 1| Trjv< (p) - "&a< (Pi) II 1 < & (8) 

By Lemma [TJ and [2j we can compute p such that J-{Tr^{p),Trj V [{pi)) = T{p,p\) and Tr j^(p) = 
Tryn(pi). Let s = i||Tr^( (p) — Trj^(pi) \\ \ and t = l\\p — p\\\. Then we have 

(R,p)-c = (R,p) + (R,p-p}-c 

1 1 

> 2" Tr - M ^ - Tr M(Pi)\\i 2 HP "Pill 



s 



t—S>s — \/2s — s 2 — S 



> 5 - yj25 -S 2 -5= -y/2S - 5 2 

where the first inequality is due to Equation QO and Lemma [5] and the second inequality comes from 
LemmaH The last inequality is because s — V2s — s 2 is decreasing when < s < 0.2 and by Equation 

□ 

Theorem 3. For any guess c (0 < c < \)for the feasibility problem, and let the input size denoted by \x\ be the size 
of function f\, there is a NC algorithm to solve the feasibility problem by returning either case in Lemma® 

Proof. The algorithm basically follows the Theorem[2]and Lemma[6] Since the particular feasibility problem 
in our consideration has Y(-) = Tr^vf (■) and A = R,B = Tr^ (pi), we have 

N(IT) = - V R + t M ®P + {pc- (Tr M ( Pl ),P))l M0V 
according to the definition in Equation (6j where EL — ^ ^ p^£T2in our problem. It is easy to see that 

||N(n)||oo < || -pR\\oo + \\1 M ®P\\oo + || (pc- (T5rjn(Pi)/P))lAf®vll<» < 1 + 1 + 1 = 3 

for any IT e T 2 and = Tr^ (p) can be calculated in NC for any p G T\. Thus, by Theorem[2j we can 

compute the approximate equilibrium value and point in Lemma[6]to precision 5 = 0(1/ polylog(\x\)) in 
NC. Based on the two cases discussed in Lemma [6[ we can either claim the original problem is infeasible or 
calculate the p by a NC algorithm according to Lemma [2] Compose all the NC circuits above, then we have 
a NC algorithm for the feasibility problem. □ 
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3.2 Solution to the Promised Version and General Case 

We are ready to apply the result of Theorem[3]to simulate QIP(2) or more general cases. Recall the definition 
of QIP(2), there will be two promises with gap A = c(|x|) — s(|x|) = Cl(l/poly(\x\)). Thus, 

Corollary 1. QIP(2) C PSPACE 

Proof. For any input x, we simply compose the following circuits. 

• For any specific x, compute the corresponding initial state p\ and the function f\. This can be done 
in NC(poly) because it only involves the computation of the product of a polynomial number of 
exponential-size matrices that corresponds to the quantum circuits used by the verifier. 

• Choose the guess value c — |(c(|x|) + s(|x|)) and precision 5 = ^ A 2 . Then use the NC algorithm 
implied by Theorem [3] to calculate the equilibrium value A* to the precision S. 

• Finally based on the two cases in Lemma [6) we can claim either the optimum value of the SDP is less 
than c or at least c — \/2S — 5 2 > c — g A. Then we are able to tell whether x £ L. 

Due to the facts in Appendix IA.21 all the circuits can be composed in NC(poly). Because of the fact 
NC(poly)=PSPACE IBor77l , we have OIP(2)CPSPACE. □ 

Furthermore, for the case where no such promise exists we can develop a binary search to approximate 
the optimum value to high precision. 

Theorem 4. Let x be any instance of SDP (I) and a be the optimum value of SDP (I). There exists a NC algorithm 
which can calculate a to precision S = Cl(l / polylog(\x\)) . Furthermore, there is a NC algorithm to compute pi 
such that (R,p2) > oc — 5. 

Proof. The proof follows from the binary search based on Lemma[6] Start with guess value c, by Theorem|3l 
we have a NC algorithm to claim either a<cora>c — S. Thus, by using binary search, we can calculate 
cc to precision S. Since S = Q(l/ polylog(\x\)), there will be at most polynomial-logarithm iterations in the 
binary search. Therefore, all the circuits above can be composed in NC. □ 

4 One-round Product QRG 

In this section, we demonstrate how the Framework [2] can be applied to real problems. Here we consider 
the simplified version of quantum refereed game with one round (two turns). The upper bound of the 
complexity class recognized by the latter model, denoted by QRG(2), becomes more and more interesting 
after the proof QIP=PSPACE |JJUW09| [M50ll. Particular, it is interesting to see whether QRG(2)=PSPACE 
while its classical counterpart RG(2) equals PSPACE IIFK97I . 

The general one-round quantum refereed game works as follows. After receiving some input x, the 
verifier then prepares some quantum messages and send them to both Yes and No provers. After both 
provers reply with quantum messages, the verifier will base on all the quantum states at his hand to decide 
whether to accept x or not. Now, let us consider a simplified case where the messages sent to the Yes prover 
and No prover are product states. We denote all the languages recognized by this procedure by product- 
QRG(2). At the first sight, this might seem to be a very restricted complexity class. However, by using the 
techniques from the recent result IBSW10|, we can prove the product-QRG(2) contains all the languages 
which can recognized by the most general model of one-round quantum refereed game except the message 
sent to Yes prover is only of poly-logarithm size. 

Let us formulate the product one-round quantum refereed game in the following way. Since the mes- 
sages sent to both provers are product states, let Vy denote the verifier's private space when interacts 
with the Yes prover and y denote the message space between the verifier and the Yes prover. Similarly, 
let Vjv denote the verifier's private space when interacts with the No prover and M denote the message 
space between the verifier and the No prover. Without lose of generality, we can assume the pure states 
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Py G D (Vy ® y) ,Pn G D (Vjv ® A/") are the initial states of the verifier given some input x. We will assume 
the input x is fixed in the following discussion. Then the Yes and No provers will apply some admissable 
quantum operations ^>y es , ^No respectively on part of the density operators py, Pn- The resultant states will 
be Cy = <5y es ® 1l(Vy)(i°y)' <7 n — ^No ® 1l(v n )(i°n)- Finally, the verifier will make some POVM measure- 
ment on uy eg) <7 N to decide whether to accept or reject. Let R (0 < R < 1) be the POVM which corresponds 
to the case where Yes prover wins, then Game Value GV(R) is defined to be 

GV(R) = max min (R, ay eg) o~n) — min max (R, <7y (g) o~n) 

where 

Ti = {cr G D ( Vy <g> y ) : 3 admissable <E> Yes : L (y) — > L (y) ,cr = <J> Yes eg l L(Vy) (py ) } 

and 

T 2 = e D(V N ® A/") : 3 admissable $ No : L (A") -> L (A/") ,(7 = 4> No <g 1 L (v n )(i°n)} 
Since py, Pn ar e pure states, by Lemma|4j we can simplify the definition of T\, T 2 to 

T\ = {<Ty 6 D (Vy eg y) : Try (cry) = Try (p Y )} , T 2 = {cr N e D (V N ® A/") : T W (cr N ) = Tr^{p N )} 

Definition 2. Any language L is inside this product-QRG(2) if and only if for any input x, 

• Ifx G L,thenGV{R) > c{\x\) 

• Ifx i L, then GV{R) < s(\x\) 

where c(\x\) — s(\x\) = £1(1/ poly (\x\)). 

Let b = j(c(\x\) + s(\x\)) and A = c(\x\) — s(\x\). By applying Framework|2 we define the convex- 
concave function h\ as follows. 

2 

%({cry,n},cr N ) = (K,c7y(g)c7 N ) - b+ - (n,Tr_^(c7 N ) -Tr^(p N )) 

where n comes from the set T — {IT : < II < ly N }- The function hi is actually the weighted sum where 
the factor a is chosen to be 2/ A . Then we will consider the new equilibrium value of function h\ instead. 
Precisely, we define 

u* = max min h\({ay ,H\ , Cn) = min max h\ ({cry, TI}, o~n) 

{(7y,n}GT 1 xTcr N eD(V N !giN) <t n gD(V n ®N) {oy,n}6Tj x T 

Then we are ready to show the gap between two promises in Definition [2] can be transferred to the equilib- 
rium value fi*. 

Theorem 5. Given the two promises in the Definition^ we have for any input x, 

• Ifx e L, then u* > ^A 

• Ifx H, then u* < -\h 

Proof. • If x 6 L, choose ({Vy, LP}, C7^) to be the equilibrium point. By Lemma [5j we have 
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By Lemma[3j there exists some ojj such that Tryy/ - (cjvf) = Txj^(p^) and ^{fi-i, (Tjj) = J r (Tr J v(tr^) / Tr_/^( ( ON)). 
Lets = 2 1| on || x andt = z ||Tr^(cr*) - Tr^(p N ) \\i , thus 

^* = (K,cr^(8>^)-6 + (R,crf®(cr^-c7^)) + -||Tr^(cr^)-Tr^( l o J v)||i 

11 1 A 2 

> 2 A " 2"^~^" 1 + a^^^n) - Tl Ar(PN)\\i = 2 " S+ A f 



A 2 ,„ n r N A 2 /„ 4 

> 2-» + 5 d-Vi r ^)a2 + s -v' 1 + ^ 
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where the first inequality comes from Lemma [5] and the second inequality is due to Lemma [3] The 



third inequality is due to the fact that mino< s <i \ — s — 1 \f \ — s 2 = \ — y 1 + ^ . The last inequality 
comes from the fact 1 — Vl + x 2 > — lx 2 for any < x < 1. 

• If x i L, by Definition |2], there exists a c"n satisfying Tr^c^) < TrA/do^) and for any cry G 7i, we 
have (R, vy <g> <7~n) — b < —j. Thus, 

^* < max /zi({c7- Y ,n},c7N) < 

{a Y ,n}es 1 xT 2 

□ 

Thus in order to tell whether x is inside L, it suffices to compute the value of /(* and make the decision 
according to Theorem [5] In the following proof, we will implicitly make use of the result in Theorem |4] 
recursively at each iteration as the solution to the oracle. 

Corollary 2. product-QRG(2)C PSPACE 

Proof. The proof of this corollary is quite similar to Corollary [T] Whenever an input x is given, we can 
calculate R in NC(poly) and approximate the equilibrium value p* . By composing all these circuits, we 
prove the whole circuit is in NC(poly) and thus in PSPACE. 

The only difference is the function h\ is not in the same form as we discussed in Theorem 12 However, 
by applying the general framework for the equilibrium value (see Appendix |A.3t , we are able to calculate 

equilibrium value }i* to high precision in NC as well. Particularly, for each c^' generated, we will choose 

cr^ to be the return value by TheoremHwhere the POVM for SDP (I) is Try N ^ N ( (lv Y g>Y ® )K(1 Vy ®y ® 

\J o~l*^ ) ) . Furthermore, we will choose nW to be the projection onto the positive eignspace of Tr_\f(cr^) — 
Tr^f(p-^). Similarly to the algorithm in Figure[IJ we will choose 

N{a Y ,n) = Tr Vy » Y ((^V ® 1v n ®n)^(V^® 1v n ®n)) + |n<g> t N - (| (n,Tw(p N )) + b)l VNm 

Since A = 0(1/ poly (\x\)), we have ||N(cy, IT) ||oo is always bounded by O(poly(\x\)). Thus the whole 
algorithm can be accomplished in NC(poly). 

□ 

5 Conclusions and Open Problems 

In this paper we demonstrate how the Framework [T] can be used to solve a class of SDPs. Moreover, by the 
examples of QIP(2) and QMAM, we demonstrate how the conversion from approximate to exact feasibility 
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can be done in our framework. The generic form in Theorem [2] also illustrates the potential of our frame- 
work to solve other SDPs. In addition, our example of product-QRG(2) illustrates how the Framework [2] 
can be used to calculate the equilibrium value of more complicated form. 

However, there are several limits and unknown facts about our two frameworks. As mentioned in |JJUW09| , 
it might be impossible to solve any SDPs in NC. Thus, we cannot hope to include all possible SDPs into 
our framework. Understanding what kind of constraints for SDPs can be solved via our framework is a 
major open problem. So far, we have positive results when the constraints are partial trace and its simple 
combination. It can be easily verified that the constraint Y = Tr^ ( U ■ If * ) for some unitary !J can also be 
solved in the similar way. 

Another open problem is whether any generalization of the multiplicative weight update method can 
be found to improve the results based on the multiplicative weight update method. The recent survey 
paper IHazlOl could provide some insights into that. 

Finally, it is still open whether QRG(2)=PSPACE while its classical counterpart RG(2)=PSPACE. 
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A Extended Preliminaries 

A.l Fundamentals of Quantum Information 

In this section, we will provide a summary of the fundamental notations and facts in quantum information. 
We assume the readers are familiar with these knowledge, and most part of this section is meant to make 
clear the terminology and well-known facts used in this paper. For those readers who are not familiar 
with these concepts, we recommend them to refer to |Bha97, KS V021 iNCOOl IWat 08 1 . Our notation basically 
follows the notation in Watrous's lecture notes I Wat08l|. 

A quantum register refers to a collection of qubits, usually represented by a complex Euclidean spaces of 
the form X = C E where E refers to some finite non-empty set of the possible states. 
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For any two complex Euclidean spaces X , y, let L (X, y) denote the space of all linear mappings(or 
operators) from X to y (L (X) short for L (X, X)). An operator A £ L (X, y ) is a linear isometry if A* A = 1% 
where A* denotes the adjoint(or conjugate transpose) of A. The set of linear isometries is denoted by 
U(X,y) then. 

An operator A £ L (X) is Hermitian, the set of which is denoted by Herm (X), if A = A* . The eigenval- 
ues of a Hermitian operator are always real. For n = dim A 7 , we write Ai(A) > A2(A) > • • • > A„(A) to 
denote the eigenvalues of A sorted from largest to smallest. An operator P £ L (X) is positive semidefinite, 
the set of which is denoted by Pos (X), if P is Hermitian and all of its eigenvalues are nonnegative, namely 
A n (P) > 0. An operator p £ Pos (A) is a density operator, the set of which is denoted by D (X), if it has 
trace equal to 1. A density operator p £ D (A) is said to be pure if it has rank equal to one. An operator 
n £ Pos (X) is a projection if n projects onto some subspace of X. Furthermore, such operators only have 
eigenvalues of or 1. 

The Hilbert-Schmidt inner product on L (X) is defined by 

(A,B) = Tr A*B 

for all A, B £ L (X). 

A super-operator(or quantum channel) is a linear mapping of the form 

Y : L (X) — > L (y) 

A super-operator Y is said to be positive if Y(X) £ Pos (30 for any choice of X £ Pos (X), and is completely 
positive if Y ® 1l(.z) is positive for any choice of a complex vector space Z. The super-operator Y is said to 
be trace-preserving if TrY(X) = Tr X for all X £ L (X). A super-operator Y is admissable if it is completely 
positive and trace-preserving. Admissable super-operators represent the discrete-time changes in quantum 
systems that, in principle, can be physically realized. 

We refer to measurements, or precisely POVM-type measurements as a collection of positive semidefinite 
operators 

{P a : a £ Z} c Pos (X) 

satisfying the constraint X^ fl£ £ P a = ^-x- Here E refers to a finite, nonempty set of measurement outcomes. If a 
quantum state represented by p £ D (X) is measured with respect to this measurement, then each outcome 
a £ £ will be observed with probability (P a ,p). 

The trace norm of an operator A £ L (X) is denoted by || A || i and defined to be 

||A||i = Tr V A* A 

When A is Hermitian, we have 

|| A ||i =max{(P -Pi,A) : P , Pi £ Pos (X) , P + Pi = l x } (9) 
The spectral norm of an operator A £ L (X) is defined to be 

|| A||oo = max{||A«|| : u £ X, \\ u \\ = 1} 
Given two positive semidefinite operators P, Q £ Pos (A"), we define the fidelity between P and Q as 

J-(P,Q) = || VPa/Q h 
When P, Q are density operators, due to Fuchs-van de Graaf inequality, we have 

l-^||P-Q||i<^(P,Q)< sji- \\\P-Q\\\ do) 

Suppose p £ D (X) is a density operator, the purification of p in A 7 ® y is any pure density operator 
mm* £ D (Af ® D/) for which Try (mm*) = p. 

The following contains the proof of lemmas shown in Section IZTl 
Proof of Lemma Q] (Please note that this lemma was originally proved in many places. The following 
proof follows the one in |JUW09| . The only reason to include this proof is because we will use it to prove 
Lemma 12) 
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Proof. First, by the monotonicity of the fidelity function under partial trace, we have for any C 2 G D(A® B) 
such that Tr B c 2 = P2 the inequality F{p\,p2) > F{ff\i0~2) always holds. Thus, it suffices to show that 
equality can be achieved. 

Let V G U (A) such that ^Jp~\ ^/piV is positive semidefinite. Since for fidelity function we have F{pi, pi) = 
II s/ply/piWli then for such a V it holds that T(p\,p2) = Tx^^/pi^/piV). Now let C = A® B and G 
A® B ® Cbe the purification of c^, in particular, | U\ } is chosen to be 

\ui) = vec(VoT) 

By rearranging the coefficients we can find aX G L (£> ® C, -4) such that vec(X) = | U\) . Since | u.\ ) is also a 
purification of pi, there must exist a linear isometry U G U (.A, £> (?) C) such that, 

X = ^plU* 

Finally, let | u 2 ) = vec(y/p~2~VU*) G A ® B ® C where V, U are obtained above respectively. It is easy to 
see that |u 2 ) is a purification of P2 in the space A® B ® C. Thus, we choose o~ 2 = Trc(|w 2 ) ( M 2|) and it will 
hold that 

^Kcr 2 ) > I (vec(^U*),vec(^p2-VU*)) \ = \ (y/plU* , y/p^VU*) \ = MVpIVM) = HP\.Pi) 

□ 

Proof of Lemma|2] 

Proof. The proof of the Lemma [1] actually gives you a way to construct such a o~2 given p\,p2,°'\- Let us 
review the important steps in the proof again with more attention to the computation of each intermediate 
quantity. 

In the first step, we need to calculate a V G U [A) such that y/px^/piV is positive semidefinite. This 
can be done by calculating the singular value decomposition of y/p\^/p2 and let V = 1 — 2P where P is the 
projection onto the subspace with negative singular values. 

The second step calculates X such that vec(X) = \u\) = vec(^/d\). This can be done by simply rear- 
ranging the coefficients in the entries of yfd\. In order to get U G U (A, B ®C), we can calculate the singular 
value decomposition of yfp{ and get the inverse( or pseudo-inverse) of ^fp{. Then U = X* {y/pi^ 1 )* ■ 

Once we have LI and V, we can easily calculate o~2 by using the formula 

a 2 = Tr c (vec(^VU*)vec(^p2-VU*y) 

Due to the fact that fundamental operations of matrix and the singular value decomposition can be done in 
NC(see Fact [1] and Fact|3]l and the fact we can compose these NC circuits easily, then we can conclude that 
o~2 can be calculated in NC given the classical representations of p\, p\ and o~\ as input. □ 

Proof of Lemma|3] 

Proof. This is only a simple application of the Fuchs-van de Graaf Inequalities(Eq lTOl ). Namely, we have 

1 - s < T{px,p2) < ~ s 2 and 1 - t < F(oi,<r 2 ) < Vl - t 2 
Given the fact F(p\,p2) = T{v\,o~2), it easily follows that 

1 - s < \/l - t 2 and 1 - t < Vl-s 2 

□ 

Proof of Lemma|4] 

Proof. This lemma is a standard fundamental fact of quantum information that follows from the unitary 
equivalence of purifications and the fact p\ is a pure state. □ 
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Proof of LemmaH] 



Proof. Given the fact about the trace norm of any Hermitian operator A in Equation[9l we have 



^lll 




n : o<n<i 



n : o<n<i 



max 2 (A, E[) 



Thus, 




For the other equation, choose A 



A, and apply the equation above, we have 



max (A'.n 

ri:0<n<i s 




Thus 



min (A, TV) 

n : o<n<i 



1 
2 



□ 



A.2 Facts on NC and parallel matrix computations 

We denote by NC the class of promise problems computed by the logarithmic-space uniform Boolean cir- 
cuits with poly-logarithmic depth. Furthermore, we denote by NC(poly) the class of promise problems 
computed by the polynomial-space uniform Boolean circuits with polynomial depth. Since it holds that 
NC(poly)= PSPACE HBor77| , thus in order to simulate the algorithm above in PSPACE, it suffices to prove 
that we can simulate the algorithm in NC(poly). 

There are a few facts about these classes which are useful in our discussion. The first fact is the functions 
in these classes compose nicely. It is clear that if / 6 NC(poly) and g £ NC, then their composition go f 
is in NC(poly), which follows from the most obvious way of composing the families of circuits. Another 
useful fact is that many computations involving matrices can be performed by NC algorithms (Please refer 
to the survey |Gat93[ which describes NC algorithms for these tasks). Especially, we will make use of the 
fact that matrix exponentials and singular value decompositions can be approximated to high precision in 
NC. We will directly cite the well-prepared form of these facts in [JJUW09] . 

Fact 1. Fundamental operations like addition, multiplication of matrices can be done in NC 

Fact 2. Matrix exponentials: there exists NC algorithms such that 

Input: An n x n matrix M, a positive rational number rj, and an integer k expressed in unary notation 



Promise: \\M\\ <k. 

Output: An n x n matrix X such that ||exp(M) — X|| < n. 
Fact 3. Singular value decompositions: there exists NC algorithms such that 
Input: An m x n matrix M and a positive rational number n. 

Output: An m x m unitary matrix U ,n x n unitary matrix V and anm x n real diagonal matrix A such 
that 



(i.e., l k ). 



M — UAV* || < n. 
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1. Initialization: Pick a fixed e < \, and let W^ 1 ) = \ x 6 L (A"), D = dim A". 

2. Repeat for each t = 1, . . . ,T: 

(a) Let the density operator pW = WW / Tr W^ f ) 

(b) Observe the loss matrix M w e L (A 7 ) which satisfies -l x < M (f) < or < M (t) < 1*, update 
the weight matrix as follows: 

t 

W( f+1 ) = exp(-£ £ M (T) ) 

T=l 



Figure 2: The Matrix Multiplicative Weights Update method. 

A.3 Multiplicative Weights Update Method 

The multiplicative weights update method introduced in Section[T]is a framework for algorithm design(or meta- 
algorithm) that works as the one shown in Fig [21 This kind of framework involves lots of technical details 
and we refer the curious reader to the survey and the PhD thesis BKal07l mentioned in the introduction. 
However, for the sake of completeness, we provide the main result which will be useful in our proof. It 
should be noticed that {MW} is the freedom we have in this framework. 

Theorem 6. Assume < < 1 for all t, after T rounds, the algorithm in Fig \2\ guarantees that, for any 

p* 6 D (X), we have 

(l- e )£(pW /M W)</p*,£Mw\ + ^ (11) 

The proof can be found in Kale's thesis |Kal07| or the appendix of IWulOI . We will then discuss how 
this method can used to solve the feasibility problem and equilibrium value following the way in Kale's 
thesis ||Kal07| with more details. 

In order to solve any feasibility problem of general form, the primal-dual method will generate a series 
of candidate solutions X^-\X^ 2 \ ■ ■ ■ , X^ T ) for T rounds. For any XW in the round t, we require an oracle 
0\ to solve the following problem 

findYWs.t. (w*{Y^)-A,X^>0,(^B,Y {t) ^<c,Y^ePos(y) (12) 

The oracle 0\ will return such a Y^ or claim such a Yw does not exist. If such a yW exists, the primal- 
dual method will generate the X^ +1 ) via the multiplicative weight update method by choosing to be a 
renormalized version of Y* (Y'*') — A. Otherwise, it stops to claim the rescaled X^' is feasible to the primal 
problem and (a, X^\ > c. If the method does not stop for T rounds, the multiplicative weight update 
method will generate an approximate dual feasible solution Y. Under certain conditions, such a Y can be 
converted into exact dual feasible solution Y such that (B, Y) < (1 + e)c. Due to the duality of SDPs, this 
implies that a < f> < (1 + e)c. 

In addition to the convertibility from approximate to exact feasibility, another difficulty in applying this 
generic method is the design of the oracle 0\ . Efficient solution to the oracle 0\ is necessary to guarantee 
an efficient algorithm for the feasibility problem. We will refer this requirement as the efficient solvability. In 
addition, we need the spectrum of Y*(YW) — A of the oracle 0\ is bounded within a small range. We will 
refer this requirement as width-boundness. 

The generic framework to calculate the equilibrium value is different in the sense of the design of the 
oracle and the use of Theorem[6] Consider the value A, 

A = minmax / (x, y) = maxmin/(x,y) 



18 



for some convex-concave function over X x Y where X, Y are convex compact sets. By convex-concave, we 
mean 

Definition 3. A function f on X x Y is convex-concave if for every y £ Y the function Vx 6 X,f y {x) = f(x,y) 
is convex on X and for every x 6 X the function My £ Y,f x (y) = f(x,y) is concave on Y. 

Consider the case when X is the set of density operators up to some factor. Again, we will generate a 
series of x^\ x^\ ■ ■ ■ , x^ £ X for T rounds. For each x^' in the round t, we require another oracle O2 to 
find approximate solution, denoted by y W, to the following optimization problem, 

max/(x( t ),y) (13) 

Then j( t+1 ) will be generated via the multiplicative weight update method. After T rounds, we can claim 
Y Yj=\ f{x^\y^) is the approximate value for A. Similarly, the oracle O2 is required to be efficient solvable 
and width-bounded. By the efficient solvability, we mean the oracle O2 can be solved efficiently. By the 
width-boundness, we mean the £00 norm of y W is bounded in some sense. 

We will conclude this section with the proof of Theorem [2] 
Proof of Theorem|2] 

Proof. First note that || N^' ||oo < r for any 1 < t < T. Thus, 

< M (f) = (N {t) +rt x )/2r < l x 

Then it is easy to see this is a typical multiplicative weights update method. Due to the Theorem|6l we 
have: 

(1 - e) t (p {T) ,M^) < (p*, t M (T A + ^ (14) 

for any density operator p* eD(X). Substitute MW = (Af( f ) + r\ x ) /2r into Equation 03] and divide both 
side by T, note that (p^\ \ < 1, then we have 

I f) ( P {T) , N< T >) <\(p\ t N (T) ) + 2re + ^ (15) 

s 



By choosing e = ^ and T 



16r 2 In D 



, we have 



1? 



According to the definition of 1 < t < T, we have 

I £ (s( P W),n«) < I £ ( S (p*),nW) + j (16) 

T=l T=l 

Choose (p*, FI*) to be the equilibrium point and substitute p* into Equation 116] , we have 

T 1 T 



I £ (s(pW),nM) < 1 £ (s(p*),nw) + s < x* + s (it) 



T = l " T=l 

where the last inequality comes from Equation HJ. 
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Since p = A Et=i P^ an d by definition of n, we have 



T=l T = l 

where the first inequality is due to the definition of TV?) and the second inequality comes from Equation 
1161 . On the other side, (p, TT) > A* by the definition of fl and Equation 0). Thus, we have 

A* < < A* + 6 

Finally we need to show that this algorithm can actually run in NC if &i = 0(1/ polylog(\x\)). Consider 
every iteration of the algorithm. The only operations involved are the fundamental operation of matrices, 
the singular value decomposition and the exponentials of matrices. Due to the Fact 1112131 they all can 
be computed in NC and the circuits for the computation can be easily composed. Since there are only 

T = 16r2 J? N iterations and r = 0(polylog(\x\)),S = O(l/polylog(\x\)), thus there will be at most poly- 
logarithm iterations respect to the input size. Therefore, the whole circuit will be NC. □ 

B QMAM case 

In this section, we will demonstrate how the Framework [JJ can be applied to the SDP of QMAM. Due to 
space limit, we will directly describe the SDP used in ]JJUW09| and denote it by SDP (II). 

SDP Problem Feasibility Problem 

maximize: (R,p) ask whether: (R,p) > c 

1 1 
subject to: Try (p) < -I4 (8> o~, subject to: Try (p) < -I4 ® a, 

p G D (A ® X ® y) ,a G D (X) p G D (A <8) X <g> y) , a G D (X) 

where K(0<R<l^)isa POVM measurement and the space A is of dimension 2. We need to design an 
algorithm to distinguish between the following two promises. Let oc be the optimum value of the SDP (II). 

Definition 4. Any language L is inside QMAM if and only if 

• IfxeL,a> c(\x\). 

• Ifx£L,a.< s(\x\). 

where c(\x\) — s(\x\) = Q(l/ poly (\x\)). 

Following the Framework [TJ we consider the feasibility problem above. Precisely, we define 

MM ,n) = (('-<M Try(pMl ^),n) as) 

where {p,a} G T a = D (A ® X <g> y) x D (X) and n G T 2 = {n : < n < 1a®x®c}- Let A^ be the 
equilibrium value of function fi, namely, 

Ao = min max h({p, o~\, Tl) = max min foilp.o-} ,IY) 

{p,cr}£T 1 neT 2 J nGr2{p <£ r}€Ti 

Base on the Theorem [TJ the value of A^ will imply whether the original problem is feasible. In order to tell 
the two promises in Definition |H we will choose the guess value c = j(c(|x|) +s(|x|)). 

Lemma 7. Given the two promises in Definition^ we have 

• IfxeL, then A£ < 0. 
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• Ifx i L, then A£ > gA 2 . 
where A = c(\x\) — s(\x\). 

Proof. • If x e L, then there exists ap e D(i®^ ^y),^ e D (A") such that (p, R) > c and Try(p) < 
O cr. This implies < 0. 

• Otherwise, let ({p*, C*}, TP) be the equilibrium point. Due to Lemma[5j we have 

A2 = M{P*> <>*}> n *) = max ( c - ( R > P*) >°} + \ II Try - 5^ ® tr|| a (19) 

By Lemma [TJ there exists a p G D (A® X ® y) such that T(\1_a ® <T,Tr y (p*)) = F(p,p*) and 
Try(p) = jl^4 ® cr. Lets = i 1 1 Try (p*) — ilx ® tr||i and t = ^ II P — P*||i- Then if t < \A, we 
have 

> c- (R,p) + (R,p-p*) +s 

> ^A-f + s 

> ^A - f + 1 - \J\ - t 2 



> l A _l A + l-\/l-iA2>iA 2 
-2 2 V 4~8 

where the first inequality is due to Equation ||T9l , the second inequality comes from Lemma|5]and the 
third inequality comes from Lemma [3] The last inequality is because t + s/l — t 2 is increasing when 
< t < \ and 1 - Vl - x 2 > \x 2 for any < x < 1. On the other side, if t > \A, by Equation f!9l, 

Af > s > 1 - Vl-t 2 > h 2 > Ia 2 

2 8 

Finally, we have A| > ^A 2 in this case. 

□ 

The only part left is to prove that we can calculate the equilibrium value An to high precision in NC. As 
the readers might notice, the set T\ is no longer a simple set of density operators but a cross product of two 
sets of density operators. However, we are still able to use a modified version of the algorithm in Figure [1] 
to solve the problem. 

Precisely, we claim the algorithm in Figure [3] will be able to calculate A| to precision 5 in NC. The proof 
is almost the same as the proof for Theorem [2] The only difference is that we need to update pW,cW 
independently and get two inequalities from each update. Then we combine them to get the final result. 
This can be done because 

/ 2 ({ P ,cr},n) = (s 1 ( |0 ),n) + (s 2 ( £ r) / n) 

where 



Let n = I r p j again, we can choose 

Ni(n) = -pR + P®ty + pct A0X0 y andN 2 (n) = -It^P 

It is easy to verify that || N\ (TV) ||co, || N^II) ||oo is bounded by 3. Thus given the gap of the equilibrium value 
Xt between two promises in Lemma[7j we can distinguish them in NC(poly) namely PSPACE. 
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1. Let e = | and T = ^L^R 

2. Repeat for each t — 1, . . ., T: 



Also let WW = l A ®x®y, = 1*. 



(a) Let = wW/TrWW,cr( f ) = vW/TrVW and let TlW be the projection onto the positive 
eigenspace of Si(p^ s> ) + S2(y( f )). 

(b) LetM^ = (Ni(nW) + rl A0X0 y)/2r, M$ = (N 2 (nW) + r\ x )/2r and update the weight 
matrix as follows: 

t 

W( f+1 ) = exp(- £ £M) T) ) 

T=l 
T=l 

3. Return j YJc=\ (Slip ) + S2{c^),Tl^ \ as the approximate equilibrium value of A|. 



Figure 3: An algorithm that computes the approximate value of A£ to precision 5. 



C Precision Issue 

The discussions on precision issue about the parallel implementations in our paper are quite similar to 
the arguments used in |JJUW09| |JUW09| |jW09l IWulOI . Again, we will make use the three facts in Ap- 
pendix |A2] and truncate the computation to sufficient precision for each step. It should be noticed that 
because of our new framework we will only use three types of operations of matrices (namely, the three 
facts in Appendix |A2|. The analysis of the precision issue in our paper then should be easier than the one 
in IJJUW09] since the latter one involves other types of operations, like the inversion of matrices. 
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